arXi v: hep-th/9811091 vl 10 Nov 1998 


Nonlinear QED Effects in Strong-Field Magnetohydrodynamics 


Jeremy S. Heyl 

Theoretical Astrophysics 130-33, California Institute of Technology, Pasadena, California 91125 

Lars Hernquist 

Lick Observatory, University of California, Santa Cruz, California 95064, USA 

We examine wave propagation and the formation of shocks in strongly magnetized plasmas by 
applying a variational technique and the method of characteristics to the coupled magnetohydro¬ 
dynamic (MHD) and quantum-electrodynamic (QED) equations of motion. In sufficiently strong 
magnetic fields such as those found near neutron stars, not only is the plasma extremely relativistic 
but the effects of QED must be included to understand processes in the magnetosphere. As Thomp¬ 
son & Blaes 0 find, the fundamentai modes in the extreme relativistic limit of MHD coupled with 
QED are two oppositely directed Alfven modes and the fast mode. QED introduces nonlinear cou¬ 
plings which affect the propagation of the fast mode such that waves traveling in the fast mode 
evolve as vacuum electromagnetic ones do in the presence of an external magnetic field [^. The 
propagation of a single Alfven mode is unaffected but QED does alter the coupling between the 
Alfven modes. 

This processes may have important consequences for the study of neutron-star magnetospheres 
especially if the typical magnetic field strength exceeds the QED critical value (Bqed ~ 4.4 x 10^^ G) 
as is suspected for soft-gamma repeaters and anomalous X-ray pulsars. 


I. INTRODUCTION 

Ultrarelativistic plasmas play an important role in energy transmission in many astrophysical settings including 
neutron-star magnetospheres, black hole accretion disks and the sources of gamma-ray bursts. Additionally, in the 
case of neutron stars, the magnetic field may exceed the QED critical value (~ 4.4 x 10^^ G) and vacuum corrections 
may affect dynamics of the electromagnetic field. 

In relativistic field theory, it is natural to study dynamics using a Lagrangian formulation. In the case of QED, 
the one-loop vacuum corrections may be summarized by an effective Lagrangian which includes both the classical 
Lagrangian and the consequences of virtual pairs If the separation of the charges comprising the plasma can 

be neglected, the dynamics of the plasma may be treated using magnetohydrodynamics. Again the determination of 
a suitable Lagrangian expedites the relativistic treatment of magnetohydrodynamics . 


II. THE ACTION 


The Lagrangian derived by Achterberg may be simplified dramatically if the inertia of the charge carriers can 
be neglected, i.e. in the extreme relativistic limit. Thompson & Blaes Q present two formulations for magnetohy¬ 
drodynamics (MHD) appropriate in this limit. Furthermore, they also examine the limit where k']_ <C eB. In this 
limit, the passing MHD wave cannot excite the charge carriers into the second Landau level, so they are effectively 
trapped along a single field line. The fermion fields are restricted to 1 -|- 1 dimensions, and they may be treated using 
the technique of bosonization. The fermion fields are replaced by a four-dimensional axion field (0) which enforces 
the MHD condition, i.e. E • B = 0. 

Thompson & Blaes Q obtain the simple action for the electromagnetic field in the presence of the relativistic 
plasma: 


S" = / d^x 
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We have the following additional definitions. 


I = and J = 


( 1 ) 
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where is the dual to the field tensor given by 
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jg completely antisymmetric Levi-Civita tensor. 
Here we will examine the modified action where 


S” = / (Tx 


AQED(d, if) + —dj 


( 4 ) 


where K = —J^. To maintain the CP and Lorentz invariance of QED, its effective Lagrangian Aqed must be a 
function of the field scalars 


i = 2(|B|2-|Ep), (5) 

if = - (4E • B)^ , (6) 


rather than of the pseudoscalar J = — 4E • B. 


III. WAVE PROPAGATION 

To study the propagation of waves through the plasma we will use the formalism of Heyl & Hernquist . We use 
the results of Thompson & Blaes Q to describe the traveling modes. Specifically we designate the external magnetic 
field by Bq and the electric and magnetic fields associated with the wave W (5E and (5B respectively. We also have 
the constraints (5E • Bq = 0 and k • (5B = 0 where k is the wave vector. Fig. ^ depicts the geometry of the propagating 
wave and defines the three Euler angles dr, 0 and (j) used to describe its configuration [^. We will take the x, y and 
z—axes to be aligned with (5B x k, k and (5B respectively. 

The definitions allow us to calculate the invariants 

I = 2[B^ + {SBf + 2Bo6B sin 0 sin 4- - {SEf] , (7) 

J = -4((5B)(,5E)cos0. (8) 


A. The Lagrange Condition 

To calculate the equations of motion of the wave we assume that the wave fields (5B and (5E and the axion field 9 
are dynamic while the external magnetic held Bq is static. Varying the action with respect to the axion held yields, 

QL 

— = 0 ^ J = -4E • B = 0. (9) 

The held 9 acts as a Lagrange multiplier to enforce the MHD condition. The equations of motion for the helds 5E 
and 5B are more complicated than in the vacuum case because here the relationship between the wave helds and 
their potential is more complicated. 


B. The four-potential 


We will designate the potential of the wave by the four-vector = [5At,5A{x,y,z)]. To within a gauge trans¬ 
formation the vector potential is given by {SA)x = —ipiyA) and {5A)y = (JA)^ = 0. 

Let us now examine the electric held, 


5E = -V{SAt) 


dSA 


( 10 ) 
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FIG. 1. The configuration of the MHD wave. N.B.: in the configuration depicted, all three Euler angles are less than zero. 
A portion of the electric field may be due to the scalar potential. We define 

(5Ei„ = 5E + ^. (11) 

at 

Let us insist that the direction of (5E is constant in time and its magnitude is a function of y and t only. Therefore, 
we have (5E = dE{y,t)x and calculate the curl of the irrotational component to obtain that (^E)^ = 0 (z.e. J = 0, 
COS0 = 0), unless we have (5E constant throughout space. We also find that 

{6E)x = + f{x,z,t). (12) 

However, since (i5E)a; and ip are functions of y and t we find that f{x,z,t) depends only on t. (i5E)j, is still uncon¬ 
strained. However we do know that (i5E)j, = R{dEi)x. We can immediately deduce the scalar potential obtaining 

SAt =-yRf{t) - R f 'tp^t{u,t)du. (13) 
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The first group of terms yields a constant background electric field in space which can vary in time. We drop this 
extra field and obtain the following four-potential, 
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where 



(17) 

The four-potential derived here is more general than that calculated by Heyl & Hernquist 

in that it contains an 


electric field such that k • (5E ^ 0 which complicates the derivation of the equations of motion. In summary we have 

((5B), = {6B)y = (5E), = 0 (18) 

and 

((5B), = V'.y, (^E), = and (5E)j, = (19) 


C. Equations of Motion 


The equations of motion are obtained by attempting to minimize the action with respect to the electromagnetic 
potential and the axion field. We found earlier that the axion field acts as a Lagrange multiplier to enforce the MHD 
condition. For the potential we obtain the following expression 


dL 




= 0 . 


( 20 ) 


Substituting the MHD action yields, 


2e^0 - 8J 


5Aq 


ED 


dK 
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ED 


dl 


pliv 


= 0 . 


We see that the axion field and QED alter the equations of motion for the electromagnetic field, 

= r. 


( 21 ) 


( 22 ) 


where 



(23) 

with 


A^2e^0 

oK 

(24) 

SAqed 
dl ■ 

(25) 

If we neglect the effects of QED on the dynamics we obtain that 


p = 2e^B • V6» 

(26) 

J = 2e^ (E X V6» - 6»,tB). 

(27) 


If we specialize to the geometry and potential described in the previous subsections and eliminate the axionic degrees 
of freedom from the equations, we find that remaining components of Eq. can be written as 
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where 


ai'^,yy + bli’,yt + Cl'i/'.tt = 0 (28) 

a3il^,yy + = 0 (29) 


where 
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|l — Q (l + R^) 


(V’.y + Bqz)R 
Boy 


[Q^l{l + R^)-1] 


Box = Bo cos 'P cos (j) 
Boy = Bq cos 'P sin (j) 
Boz = Bq sin 0 sin 'h 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
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D. The Classical Limit 

If we neglect the effects of QED in the equations above, we can derive the wave equations in relativistic MHD. 
Specifically, we take Q = 0 (in the classical limit cIAqed/ 9/ = —1/4) to obtain 
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IM 
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(40) 

(V'.i^.yt - ^.y^.tt - Bozi’,tt) = 0 

(41) 


where 


= (42) 

Boy 

Eq. M is satisfied by 

tp{y,t) =tp{y±vt). (43) 

This equation also satisfies Eq. ^if i? = 0 or Bqz = 0; consequently, MHD supports fully nonlinear modes if dE-k = 0 
(fast modes) or dB • Bq = 0 (Alfven modes). 

For the fast modes, we have R = 0 and a dispersion relation of = fc^. For the Alfven modes, we find that Bq, dE 
and k all lie in the plane perpendicular to dB. Since dE • Bq = 0 we have 


(^E), 


Box ,2 _ 1 

~ l + i?2 


= sm 


( 44 ) 


yielding a dispersion relation of uj = ksg is the component of k directed along the external magnetic field. 

However, even if dB • Bq = 0 unless Bq and k are parallel {i.e. R — 0), 
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( 45 ) 


^(2/: t) =1pl{y + VAt) + ^’ 2 ( 2 / - VAt) 

does not satisfy Eq. 1^. Two oppositely traveling Alfven modes will interact through the cross term, 


Rv^ 

^ {i^l,Vy'lp2,y + V'2.yy'01.y) = 0. 




Oy 


(46) 


It is also straightforward to derive equations Eq. E and Eq. E from Eq. ^ and Eq. E using Maxwell’s equations. 


E. The Effects of QED 


QED introduces several additional terms into the equations of motion. By restricting the geometry of the wave, we 
can still satisfy Eq. ^ identically. 


1. Fast Modes 


First, we examine the fast modes which have both R = 0 and Bq^ = 
effects are included. For these modes, we obtain 


0 and therefore satisfy Eq. 


even when QED 


6i = (V'.y + Bo.) (47) 

Cl = 1 - (48) 

Because Oi does not depend on R or Bq^ it is still given by Eq. As in Heyl & Hernquist [H we expand the 
coefficients to first order in the fields. 


ai = ai.o + + ai.Et/’.t + 0{dB^) 

h = 5i,o + + 0{6B^) 

Cl = Cqo + Ci^Bip,y + CpB'i/’,* + 0{5B'^). 

Notice the sign change relative to Heyl & Hernquist |^. In the previous work, we selected a gauge where {5E) 
We obtain 


and 


01,0 = —4 

ai,B = -16 

bi.E = 32 
ci,o = 4 


,d2 5^Aqed 

QJ2 4 

S^Aqed ^ 

o——doQz 


dP 


52 A, 


QED 


dP 

9Aqed 




Oz 


dl 

T ^ 5^Aqed ^ 
ci,B = 16 Bo., 


9Aq 


ED 


dl 

I 9^Aqed 

dP 


Bl. 


(49) 

(50) 

(51) 

-V’.t- 

(52) 

(53) 

(54) 

(55) 

(56) 


ai,E = 5i,o = 61 ,B = ci,B = 0. (57) 

If we substitute these expansions back into Eq. we obtain an identical result to that of Heyl & Hernquist in the 
limit of AT = 0. This is not particularly surprising because the fast modes do not carry current and cannot excite the 
axion field. Furthermore, the results of Heyl & Hernquist ^ may be applied directly to understand the evolution of 
the fast modes including the effects of QED to one-loop order. 

Heyl & Hernquist Q found that a electromagnetic wave traveling through a strongly magnetized vacuum will 
develop discontinuities after traveling a distance proportional to its wavelength and amplitude. Furthermore, the 
opacity to shocking peaks near the critical field (~ 4.4 x 10^^ G). After the discontinuity forms, the energy of the 
wave is quickly dissipated most likely as electron-positron pairs. The results of this section show that the presence of 
the plasma does not affect the development of this nonlinearity for the fast modes; therefore, even in the plasma-filled 
magnetosphere surrounding a neutron star, one would expect shocks to develop as waves in the fast mode propagate. 
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2. Alfven Modes 


The treatment of the current-carrying modes within QED is more subtle. We can use the results that Bqz = 0 and 
R = —Box/Boy to simplify the equations, 


Ol = -1 - Qi^'^y 

h = (1 + R^) 

Cl = — Q^'^t (l “f R?^ -f 1 -f R? 

W. - Bl) 

Boy 

idOy 

To proceed we combine Eq. and Eq. ^ into a single relation by eliminating 4’,yt, 


where 


- -f V',tt = 0 

aibs - ash 
63 C 1 - 61 C 3 

1 1 - QV'^y 

1 + i?2 1 - Qi/;2 (1 + i?2) ■ 


(58) 

(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 


This elimination is impossible in pure MHD {i.e. Q = 0) since the term ijj^yt appears only in the second of the equations 
of motion (Eq. I)- 

Let us attempt to use a simple wave solution similar to that of MHD, 


If we substitute this into Eq. |6^ , we find 

^A,QED 


= CA,QEDi) ■ 


1 _ 1 - _ 

l + R^l- Qci (1 -f i?2) 


(67) 


( 68 ) 


which is satisfied by 


^A,QED - ^'A - 


(69) 


We find that a single Alfven mode does not suffer any nonlinearities due to QED if the QED effective Lagrangian is 
invariant with respect to gauge and Lorentz transformations and 


^ ^ 0 . (70) 

Eq. ^ like Eq. indicates that two oppositely traveling Alfven modes will interact. If we substitute Eq. ^ into the 
equations of motion for Alfven waves including QED we obtain 

Qi’ly^l.yy + {1 ^ 2} = 0 (71) 

73,, 2 

-j^Mi,yy {l + Q[Bl + 2^2,y (V'l.y - ^2,y)] } + {1 ^ 2 } = 0 (72) 

OOy 


where {1 2} designates the same terms repeated with the functions ifji ^nd ' 4’2 swapped to obtain a symmetric 

sum. 

To lowest order in the wave helds we find that QED does not affect the coupling between the Alfven waves. However, 
at higher order, vacuum processes introduce several new interaction terms. 
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IV. CONCLUSIONS 


We find that QED affects the propagation of magnetohydrodynamic fast modes traveling through an external 
magnetic field. The induced nonlinearity is identical to that suffered by electromagnetic radiation traveling through 
an external magnetic in the absence of a plasma. For the Alfven modes, we find that QED introduces additional 
couplings between oppositely directed waves, but that because of the gauge and Lorentz invariance of QED, no 
nonlinearities manifest themselves in the propagation of a single Alfven mode. 

These new nonlinear processes emerge in regions where the magnetic field strength is comparable to or exceeds 
Bqeb ~ 4.4 X 10^^ G. Both anomalous X-ray pulsars (AXPs) and soft-gamma repeaters (SGRs) circumstantially 
exhibit such strong magnetic fields consequently, these nonlinear processes specific to intense magnetics fields 

may play a important role in the magnetospheres surrounding these objects. 
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